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STABILITY  OP  FLOW  BETWEEN  ARBITRARILY  SPACED  CONCENTRIC  CYLINDRICAL 
SURFACES  INCLUDING  THE  EFFECT  OF  A  RADIAL  TEMPERATURE  GRADIENT 

by 

* 

J.  Walowit,  S.  Tsao,  R.  C.  DiPrIm 


ABSTRACT 

The  stability  of  Couette  flow  and  flow  due  to  an  azimuthal  pressure 
gradient  between  arbitrarily  spaced  concentric  cylindrical  surfaces  is  in¬ 
vestigated.  The  stability  problems  are  solved  by  using  the  Oalerkin  method  in 
conjunction  with  a  simple  set  of  polynomial  expansion  functions.  Results  are 
given  for  a  wide  range  of  spacings.  For  Couette  flow,  in  the  case  that  the 
cylinders  rotate  in  the  same  direction,  a  simple  formula  for  predicting  the 
critical  speed  is  derived.  The  effect  of  a  radial  tenperature  gradient  an  the 
the  stability  of  Couette  flow  is  also  considered.  It  is  found  that  positive 
and  negative  tenperature  gradients  are  destabilising  and  stabilising  respec¬ 
tively. 


INTBODUCTIDN 


The  stability  of  viscous  flow  between  concentric  rotating  cylinders 
was  first  investigated  experimentally  and  theoretically  by  G.  I.  Taylor  in  his 
famous  paper  in  1923  [1].*  He  showed,  in  agreement  with  experiment,  that  when 
a  certain  criterion  was  violated  a  secondary  motion  in  the  form  of  toroidal 
vortices  sets  in.  For  the  most  part  Taylor's  computations  were  limited  to  the 
case  when  >0  or  only  slightly  negative  (with  a  few  exceptions) 

and  Vj  “  R3/R2  ~1  where  and  n2  are  the  radii  and  angular  velocities 

of  the  inner  and  outer  cylinders  respectively.  However,  he  did  give  a  correc¬ 
tion  formula  suitable  for  1-Vj  small.  Since  this  original  investigation, 
there  have  been  a  number  of  papers  dealing  with  suitable  mathematical  techniques 
for  solving  the  related  eigenvalue  problem  for  the  "small-gap"  problem  ( VJ  — ► 1 ) , 
but  for  arbitrary  rotation  rates.**  The  case  of  -Jl  large  is  particularly 
difficult.  (See  [28],  [2]).  On  the  other  hand,  the  "finite-gap"  problem 
(  /  1)  has  only  recently  been  considered.  One  of  the  difficulties  in  treating 

the  finite-gap  problem  is  that  the  differential  operators  in  the  eigenvalue 
problem  has  variable  coefficients  in  contrast  to  the  constant  coefficient  opera¬ 
tors  that  appear  in  the  small -gap  problem. 

Using  an  expansion  technique  similar  to  that  for  the  small-gap  problem 
in  which  one  of  two  equations  iB  solved  exactly,  Chandrasekhar  [3]  has  consi¬ 
dered  the  finite-gap  problem.  The  computations  are  rather  tedious  and  numerical 


*Humbers  in  brackets  designate  References  at  the  end  of  the  paper. 

**For  a  discussion  of  some  of  this  work  as  well  as  references  to  other  papers, 
the  reader  is  referred  to  Chandrasekhar  [2,  Sec.  71]. 


results  were  obtained  only  for  the  case  Xj  -  i,  but  for  a  wide  range  of  values 
of  yU,  -0.5  4  0.25.  (The  flow  is  stable  for  JU.  >  fj2".)  More  recently 

Chandrasekhar  [U]  and  Chandrasekhar  and  Elbert  [53  have  shown  that  sane  of  the 
numerical  difficulties  can  be  overcome  by  considering  the  corresponding  adjoint 
eigenvalue  problem.  Witting  [6],  treating  1-rj  as  a  small  parameter,  has  also 
considered  the  stability  problem  by  expanding  in  powers  of  (l-/j  ).  Finally 
Kirchgassner  [7]  has  constructed  the  Greens  function  for  the  finite-gap  problem, 
and  has  solved  the  resultant  integral  equation  by  an  iteration  technique.  Re¬ 
sults  were  obtained  for  i  £  VJ  <  1  with  yU*  0,  -0.U  ^  0.25  far  Fj  ■ 

and  for  -O.U  V9  for  Yj  -  2/3. 

A  similar  type  of  instability  problem  occurs  when  a  viscous  fluid 
flows  in  a  curved  channel  under  a  pressure  gradient  acting  round  the  channel. 

The  small  gap  form  of  this  stability  problem  was  first  considered  by  Dean  [83. 
The  finite-gap  Dean  problem  has  not  been  considered  in  the  literature.  How¬ 
ever,  Meister  [93  has  considered  the  more  difficult  finite-gap  combined  Taylor- 
Dean  problem  for  the  cases  fj  -  3A  and  5/6.  By  properly  rescaling  his  results 
it  is  possible  to  obtain  approximately  the  criterion  for  the  Dean  stability 
problem  for  these  values  of  /| .  However,  there  is  Insufficient  data  to  de¬ 
termine  this  limit  accurately. 

In  the  present  paper  the  finite-gap  Taylor  and  finite-gap  Dean 
problems  are  considered.  It  is  shown  that  each  of  these  problems  can  be 
solved  using  the  Qalerkin  method,  the  eigenfunctions  being  represented  as  a 
combination  of  simple  polynomials.  With  this  choice  of  expansion  functions 
all  necessary  integrals  can  be  easily  evaluated  and  numerical  computations 
have  been  carried  out  for  a  much  wider  range  of  rj  (0.1  4  xj  4  l)  than  in  previous 
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computations .  The  results  are  in  good  agreement  with  the  known  theoretical 
results  and  all  of  the  available  experimental  results. 

Finally  the  effect  of  a  radial  temperature  on  the  stability  of  Couette 
flow  between  rotating  concentric  cylinders  is  considered.  The  corresponding 
theoretical  small-gap  problem  has  been  considered  by  Yih  [lO] ,  Lai  [ll] ,  and 
Becker  and  Kaye  [12].  Experimental  results  for  gases  bave  been  obtained  by 
Becker  and  Kaye  [13]  and  Bjorklund  and  Kaye  [20],  and  for  liquids  by  Nissan 
and  Haas  [29]  and  Nissan,  Ho  and  Nardacci  [30].  With  just  a  few  additional 
computations  this  problem  can  be  solved  by  the  Galerkln  method  used  in  solving 
the  Taylor  and  Bean  stability  problems. 

THE  STABILITY  PHDHLEH3 

Let  r,G,z  denote  the  usual  cylindrical  coordinates,  with  the  %  axis 
coinciding  with  the  axis  of  the  cylinders;  and  let  RpRg  and  XI  ^  and  XI 2  de¬ 
note  the  radii  and  angular  velocities  of  the  inner  and  outer  cylinders  respec¬ 
tively.  If  ur,  Uq,  and  uz  denote  the  components  of  velocity  in  the  increasing 
r,  6,  and  z  directions  and  p  denotes  the  pressure,  the  Navi er -Stokes  equations 
admit  a  steady  solution  of  the  form 

«r  “  "  0  ,  ug  -  V(r)  ,  -iWbr  -  $V2/r  (1) 

where  f  is  the  density. 

Now  superimpose*  on  this  steady  motion  a  small  rotationally  sym¬ 
metric  disturbance  of  a  form  such  that  the  6  component  of  the  velocity  is 

Ug(r,z,t)  -  V(r)  +  v(r)  e3"*1  cosXz  .  (2) 

*A  complete  derivation  of  the  disturbance  equations  can  be  found  in  Chapter  7 
of  reference  [2]. 
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The  motion  will  be  unstable  if  there  exist  solutions  of  the  resultant  eigen¬ 
value  problem  for  which  the  real  part  of  T  is  greater  than  seroj  it  will  be 
stable  if  all  the  solutions  have  real  part  of  <T  less  than  eero.  It  is  known 
from  experimental  evidence  that  not  only  is  the  Instability  of  this  spatial 
form,  but  that  also  a  new  steady  secondary  motion  occurs.  Hence  we  can  an¬ 
ticipate  that  the  marginal  state  is  given  by  <T  ■  0  rather  than  just  real 
part  <T  -  0.  In  this  case  the  linearised  equations  for  the  distarbanoe  veloci¬ 
ties  are 


(DD*  -  X2)2  u  •  -ft  (r)  v 

(3) 

(DD*  -  X2)  ▼  “  -i-(D*V)  u 

(U) 

where  u  is  the  disturbance  velocity  in  the  radial  direction, 

matte  viscosity,  /l(r)  -  V(r)/r,  and 

i>  is  the  kine- 

D  -  4-  ,  D*  -  D  +  i 

dr  *  r 

(5) 

The  requirement  of  no  slip  at  the  boundaries  gives  the  boundary  conditions 

u  ■  v  ■  Du  •  0  (6) 

at  r  ■  and  r  -  We  will  consider  two  examples. 

The  Taylor  stability  problem.  In  this  case  the  basic  motion  is  that 
due  to  the  motion  of  the  cylindrical  walls  (Couette  flow),  and  ?(r)  is  given 
*7 

V(r)  -  Ar  ♦  | 
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(7) 


where 


Art  - 


A  -% 


■t-rt 


(8) 


lotice  that  D*V  ■  2A,  a  constant. 

It  is  convenient  to  choose  our  length  scale  so  that  the  range  of  the 
independent  variable  is  independent  of  the  parameters  of  the  problem.  Further, 
the  natural  scale ,  except  in  the  limiting  case  JA  ■  2 / (See  C28] » ) , 

for  \  is  the  gap  width.  Hence  we  introduce  the  dimensionless  variables 

r  ■  RQ  +  dx  ,  a*  X.d  (9) 

where  d-Rj-R-Landl^  -  (Rg+R^/Z.  Then  the  Taylor  stability  problem  takes 
the  form 

(DD*  -  a^)2  u  *  -  a^  T  g(x)  V  (10) 

(ED*  -  a^)  v  -  u  (11) 


where 

il(r)  -  .Q^  g(x)  ,  g(x)  -  ^ 

)miqm  V  (1- '))(***>  »  >]-V*2  »  /- ivni  (12) 


and  T  is  the  usual  Taylor  number 

T  -  -  UiXl-d^/p2  (13) 

1  r  fit 
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I 


and  we  now  use  D  and  D*  for  the  dimensionless  operator* 


”4 


5  J 


Finally  u  has  been  replaced  by  (2Ad2/i>  )“^  u* 

Equations  (10)  and  (11)  with  the  boundary  oonditii 


u  -  v  -  Du  -  0  (15) 

at  x  »  deteradne  an  eigenvalue  problem  of  the  font 

F(T,a,yU,»j)  -  0  .  (16) 

For  fixed  values  of  JJL  and  y  ,  which  determine  the  geometry  and  basic  velocity 
distribution  up  to  scale  factors,  we  wish  to  determine  the  minimum  positive  T 
over  all  real  positive  a  for  which  there  exists  a  non- trivial  solution  of  the 
eigenvalue  problem.  The  critical  value  of  T  determines  the  critical  speed  D~  ^ 
for  fixed  JX.  and  rj  ,  and  the  corresponding  critical  value  of  a  determines  the 
spacing  of  the  vortices  in  the  z  direction. 

It  should  be  noted  that  the  Taylor  ..umber  T  as  defined  by  Iq.  (13) 
is  different  from  that  used  by  Chandrasekhar  [33  (T  -  -  liABR|/p  2)  who  used 
R2  as  his  reference  length.  The  present  choice  is  the  usual  Taylor  number 
for  the  small-gap  problem;  hence,  in  the  small-gap  limit  rj-wl,  Eqs.  (10) 
and  (11)  reduce  to  the  usual  small-gap  equations.  Further  we  note  that  in 
the  case  JU  ■  0,  which  is  of  most  engineering  Interest,  both  Arj,  and  are 
proportional  to  2.  Hence  to  obtain  a  meaningful  problem  for  Jj.  ■  0  as  *  0 
it  is  necessary  to  rescale  T  by  a  factor  of  ^  c  treating  T  3 c  as  the  new  eigen¬ 
value.  That  is,  we  can  anticipate  that  for  JX.  -  0,  T-»C/yj  2  where  C  is  a 
constant  as  "rj  -»0.  Most  of  our  results  are  given  in  terms  of  T rj 2  which  also 
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reduces  to  the  snall-gap  Taylor  number  as  rj  -+1. 

The  Dean  stability  problsa.  In  this  case  the  cylinder  nails  are  at 
rest  and  the  notion  is  due  to  a  constant  cireunferaatial  pressure  gradient. 

The  basic  Telocity  is 

T(r)  -  -TT  4*(r  la  r  ♦  Or  ♦  S)  (17) 

2  fd  «8  r 


where 


R|  In  R2  -  B£  In 

4-4 


,  R2 
In  — 


■4-"?  ^ 


With  z  and  a  defined  by  Sqs.  (9),  and  T  ,  D  and  D*  defined  by  Iqs. 
(lU),  the  finite-gap  Dean  stability  problan  takes  the  foxn 


(DD*  -  a2)2  u  -  a2  p[h(x)/5]  ▼ 
(ED*  -  a2)  ▼  -  (D*h)  u 


where 


*2  ?P  (1-n2)2  -  U>)2(ln>l  )2 


7(r)  -  7  h(x)  ,  7  -  . 

•  "  2J*>  >•  Ud-rjHi-f) 


*)  * -  2  1  < - — j  T  Ul-  h2)  lnt  ♦  *2  In  >1  (1  -  1  ) 

l*f(l*  )“  -  (1-  >  l  I  )  I  !  T 


Vd  9  j 

P  “  2(“J-)  ^ 


and  7^  is  the  nean  Telocity.  Finally  u  has  been  replaeed  by  (7,^/j  )**1  u. 
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The  boundary  conditions  are  given  by  Eqs.  (15). 

In  the  small-gap  limit,  Yj  1,  Eqs.  (19)  end  (20)  reduce  to  the 
usual  small-gap  equations  within  an  integral  multiplier  of  P. 

METHOD  OP  SOLUTION 

The  eigenvalue  problems  for  the  finite-gap  Taylor  and  Dean  insta¬ 
bility  problems  can  be  conveniently  solved  by  the  Qalexfcin  method.  The  func¬ 
tions  u  and  v  are  expanded  in  complete  sets  of  functions  satisfying  the 
boundary  conditions  (15).  The  coefficients  in  the  series  are  determined  by 
the  requirement  that  the  error  in  Eqs.  (10)  and  (11),  for  the  Taylor  problem, 
or  in  Eqs.  (19)  and  (20),  for  the  Dean  problai,  be  orthogonal  to  the  expansion 
functions  for  u  and  v  respectively.  This  leads  to  a  system  of  infinitely  many 
linear,  homogeneous  equations  for  the  coefficients  in  the  series.  For  a  non¬ 
trivial  solution  it  is  necessary  that  the  determinant  of  the  system  of  equa¬ 
tions  vanish,  and  this  gives  a  detemlnantal  equation  for  T(ayU,>j  )  or  P(a,  rj  ). 
In  practice  only  a  finite  number  of  terms  are  used  in  the  series  for  u  and  v, 
say  M,'  and  this  leads  to  a  determinant  of  size  2M.* 

In  principle  there  are  a  number  of  possible  sets  of  complete  functions 
which  can  be  used  for  the  expansion  of  u  and  v.  For  example,  the  set  of  func¬ 
tions  sin(2n-l)7|’x  and  cos2nTpc  (n-1,2,. .. )  could  be  used  for  v(x),  and  the 
set  of  functions  Cn(x)  and  Sjj(x)  tabulated  by  Harris  and  Reid  [323  could  be 
used  for  u(x).  Unfortunately,  the  choice  of  tri gnome trie  functions  leads  to  a 
prohibitive  amount  of  labor  just  to  evaluate  the  entries  that  appear  in  the 
determinantal  equation.  Indeed  it  is  not  clear  that  many  of  the  integrals  can 

# 

It  is  not  actually  necessary  that  the  number  of  terma 
v  be  the  same. 


in  the  series  for  u  and 


•ran  be  evaluated  In  closed  form. 

For  our  purposes  it  is  more  convenient,  and  equally  as  satisfactory, 
to  use  simple  polynomials  in  x  that  satisfy  the  boundary  conditions.  Thus  we 
choose 

90 

u(x)  -  ZL  a<n',n(x)  »  ▼(*)  "  ZL  6  nTn(*)  (22) 

n-1  n-1 

A 

where  %(*)  “  (x2  -  J)  x11-1  and  vn(x)  -  (x2  -  a?1*1  .  These  polynomial 

expansion  functions  have  also  been  used  by  Kurzweg  [3l]  in  his  treatment  of 
several  small-gap  problems  in  hydrodynamic  and  hydroma gnetic  stability.  It 
is  particularly  interesting  to  note  that  with  this  choice  of  expansion  func¬ 
tions  all  of  the  various  integrals  can  be  expressed  exactly  as  polynomials  in 
or  as  the  products  of  In  rj  and  polynomials  in  rj  .  In  the  use  of  the 
Qalerkin  method  for  these  problems  it  is  helpful  to  note  that  the  operator 
d/dr(d/dr  ♦  r"1)  is  a  Stum-Liouville  type  operator  with  a  corresponding  weight 
function  r,  or  in  our  dimensionless  variables  "J  .  Thus  we  require,  for  example 
in  the  Taylor  problem,  that  the  error  in  the  equation  for  u  be  orthogonal  to 
^u^x),  m-l,2,...,M,  which  gives 

{*n(K»<DD*  "  #2)2  0  +  P n^O  V^n)}  “  0  » 

m  -  1,2,....  (23) 

where  (%,%)  denotes  the  integral  of  u,,,^  from  x  ■  -  £  to  x  -  J.  Similarly 
for  the  equation  for  v(x)  and  for  the  Dean  stability  problem.  The  number  of 
computations  that  are  required  are  considerably  reduced  by  the  syaaetries 
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( ^%,(DD*-a2)2  %)  -  ( )  (DD*-a2)um),  etc.  and  by  the  fact  that  the  differ¬ 

ential  operators  for  the  Taylor  and  Dean  stability  problems  are  the  same. 

As  mentioned  earlier,  while  it  is  possible  to  evaluate  all  integrals 
exactly,  it  is  actually  more  convenient  in  practice  to  proceed  in  a  slightly 
different  manner.  The  details  of  the  evaluation  of  the  integrals  are  given 
in  the  Appendix.  The  determinant  is  evaluated  by  the  method  of  pivotal  con¬ 
densation,  the  zero  being  found  by  trial  and  error.  All  of  the  computations 
were  carried  out  on  an  I. B.M.  6$0. 

For  the  Taylor  stability  problem  computations  were  carried  out  for 

M  *  2,  3,  and  It  for  a  range  of  values  of  Y\  from  1  to  0.1  and  for  a  range  of 
2  2 

jk-  <  •  For  JX.  >  the  flow  is  known  to  be  stable  to  disturbances  of  the 


fj  *  1,  and  JX 


in  the 


form  (2).  The  results  are  tabulated  in  Table  1.  For 
range  -1  £ JXftl  the  results  for  M  ■  U  are  within  0. 2%  of  the  "exact  results" 
given  by  Chandrasekhar  [?,  Sec.  7l] .  Indeed  even  with  M  »  2  the  error  for  the 
worst  case,  Jk  *  -1,  is  only  $%.  If  we  bear  in  mind  that  the  analytical  and 
numerical  work  using  the  Qalerkin  method  with  M  -  2  is  much  less  than  that  for 
the  expansion  procedure  of  Chandrasekhar  or  the  integral  equation  technique, 
this  result  is  surprisingly  good.  For  rj  •  $  and  a  range  of  jk,  -  %  ijki 
and  M  ■  1*  the  maximum  difference  between  the  present  results  and  those  of 
Chandrasekhar  and  Elbert  occurs  at  fj  -  -  and  is  about  2%.  With  de¬ 
creasing  /j  (wider  gap)  and  decreasing  jk.  it  becomes  more  difficult  to  approxi¬ 
mate  the  eigenfunction,  so  that  at  VJ  m  0.2  it  was  not  possible  to  determine  T 
accurately  for  negative  values  of  JX.  without  taking  more  terms  in  the  series 
for  u  and  v.  By  comparing  the  results  for  M  ■  2,  3,  and  U  it  is  estimated  that 
all  the  results  given  in  Table  1  for  M  -  1*  are  correct  within  2%.  For Jk.  -  0, 


11 


and  a  range  of  Yj  ,  0.2  <  YJ  <  1,  soae  recant  results  by  Roberts  [lkl  using 
direct  numerical  integration  are  essentially  Identical  with  those  given  here. 

For  the  Dean  stability  problem  computations  have  been  carried  out 
for  a  range  of  Vj ,  0.1  $  >j  $  1.  The  results  are  tabulated  in  Table  2.  The 
only  known  results  for  this  problem  are  for  Vj  ■  l.  in  this  case  the  result 
for  (VKd/\>  )  (d/Rj)^  is  essentially  identical  with  the  value  35.9k  reported 
by  Reid  [l5]  and  the  values  given  by  HSmmerlin  Ll6l»  and  Dean  [83.  For  M  »  2 
the  difference  is  only  2.6%.  Comparison  of  the  results  for  M  »  2  and  k  indi¬ 
cates  that  the  error  increases  with  decreasing  Yj,  the  percentage  change  being 
about  T%  at  Vj  -  0.1.  It  should  be  noted  that  the  percentage  change  in  the 
actual  eigenvalue  P  appearing  in  Eq.  (19)  is  twice  that  in  (VBd/ 1)  )(d/Rg)^. 

EFFECT  OF  A  RADIAL  TEMPERATURE  GRADIENT 

In  this  section  we  will  consider  the  effect  of  a  radial  temperature 
gradient  on  the  finite- gap  Taylor  stability  problem.  The  presence  of  the  radial 
temperature  gradient  will  give  rise  to  convective  effects  througi  the  interac¬ 
tion  of  the  associated  radial  density  gradient  with  the  radial  acceleration. 

In  this  analysis  we  will  neglect  the  effect  of  gravity.  Since  in  both  gases 
and  liquids  the  centrifugal  force  field  vfyr  corresponding  to  the  critical 
speed  is  often  much  less  than  g,  this  is  a  serious  approximation.*  However, 
for  stationary  cylinders  no  significant  changes  in  the  heat  transfer  as  pre¬ 
dicted  by  conduction  were  observed  in  the  measurements  of  Krausshold  [rQ  for 
Rayleigh  numbers  (based  on  g)  less  than  1(P.  For  horiiontal  cylinders  this 
was  verified  by  the  experiments  of  Liu,  Mueller,  and  Landes  [l8j  and  the 


*The  ratio  (B^ fl^/R^g  is  given  approximately  by  T(\)2/d3)jj(l+Vj)Ag(tj2-yU) 
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theoretical  analysis  of  Crawford  and  Lenlich  [19] .  In  addition,  in  the  experi- 

4 

mental  stability  work  of  Becker  and  Kaye  [13]  and  Bjorklund  and  Kays  Ij2Q]  with 
gases  no  significant  effect  of  gravity  was  observed.  Further  we  mi^it  note 
that  in  stability  problems  of  this  type  the  critical  Taylor  ntmber  is  rather 
insensitive  to  variations  in  the  basic  velocity  and  teiqperature  profiles.  Thus 
we  will  ignore  the  apparently  weak  convective  motions  due  to  gravity  in  the 
present  analysis.  The  effect  of  essentially  a  radial  gravity  field  on  various 
curved  flows  has  been  considered  by  GBrtler  [2l]  and  Kirchgtfssner  [22]. 

In  addition  to  the  three  momentum  equations  and  the  continuity  equa¬ 
tion  we  must  now  include  an  energy  equation  and  an  equation  of  state  in  our 
analysis.  We  shall  assume  that  our  fluid  has  constant  thermal  conductivity, 
specific  heat,  and  viscosity*,  and  we  will  neglect  viscous  dissipation  and 
energy  associated  with  change  in  pressure  in  the  energy  equation.  With  these 
assumptions  the  basic  azimuthal  velocity  distribution  is  given  by  Iq.  (7),  and 
the  temperature  distribution  (*)  is  given  by 

©<r>  ■  I,  -  (Tj-Tj)  (2M 

where  T^  and  T2  are  the  temperatures  of  the  inner  and  outer  cylinders  respec¬ 
tively.  As  an  equation  of  state  we  assume  ^  -  ^0[l-o<(  ®  -  ®Q)1  where  (h)  q 
is  a  reference  temperature  and  o(.  is  the  bulk  coefficient  of  thermal  expansion. 

We  now  consider  rotationally  symmetric  disturbances  of  this  steady 


*Thus  our  analysis  is  primarily  limited  to  gases  with  small  temperature  dif¬ 
ferences  or  liquids  with  only  slight  viscosity  variation  with  temperature, 
the  corresponding  small  gap  problem  for  liquids  with  variable  viscosity  has 
been  considered  by  Walowit  £23]. 
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▼elooity  and  temperature  distribution.  A  snail  disturbance  6(r,z,t)  in  the 
temperature  will  lead  to  a  disturbance  in  the  density.  For  small  temperature 
variations,  the  variation  in  density  will  be  snail  and  we  will  assume  that 
the  density  is  a  constant  ^  ,  except  when  multiplied  by  the  centrifugal  ac¬ 
celeration  term  V2/r.  Here  the  density  disturbance  is  given  by  -  ^o(©' .  This 
is  the  classic  Boussinesq  approximation.  The  disturbances  can  be  expressed  in 
the  form  given  in  Eq.  (2),  that  is,  0»  ■  ©(r)  e^*  cosX*  and  we  obtain,  using 
the  principle  of  exchange  of  stabilities,  the  following  eighth  order  system  of 
equations 

i)(ED*  -  X2)2  u  -  X2(2A*  ej 

0(DD*  -  X2)  ▼  ”  (2A)  «  (2$) 

K(D«D  -  X2)  e  *  (D®)  u 

with  the  boundary  conditions 

U“Du“V“9"0  (26) 

at  r  *  %  and  R2.  Here  K  is  the  thermal  diffueivity.  In  dimensionless  varia¬ 
bles  Iqs.  (2$)  take  the  form 

(DD*  -  a2)2  u  ■  -  a2T  jg(x)  v  ♦  N  J  [.g(x)J2  ©] 

(DD*  -  a2)  v  ■  u  (27) 

(D*D  -  a2)  6  *  u/j 

where  x,a,g(x), J  ,T,D,  and  D*  are  defined  by  Iqs.  (9),  (12),  (13)  and  (lU). 
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Also  u  has  btm  replaced  by  (2Ad2/i>  )-1  u  and  0  has  been  replaced  by 
[pr(T2-|r3i)/2iEf  Jn(l/yj  )]  8,  where  Pr  -  /  K  is  the  Fraadtl 

Pr  •(.(Tg-T,/ 

N  -  —5 - 

•*L 


is  the  ratio  of  a  Rayleigh  number,  Ray  -  C<(T2-Ty)d3  Rj rnlth 
Hyj  -  -  d/ln*j ,  based  on  the  centrifugal  acceleration  «*d  the  Taylor 

nunber  T. 

In  the  small-gap  limit,  yj  -ml,  N-*-Pr<X(T2-T-L)/2(  1-^/0.  Moreover, 
M  y|-^i,  ^  -ml,  D*“ mD;  hence  it  follows  from  the  last  of  Eqs.  (27)  and  the 
boundary  conditions  (26)  that  6  -  v.  Consequently  the  eighth  order  system  of 
equations  reduces  to  a  sixth  order  system.  This  problem  has  been  considered 
by  Tih  [jo]  and  Becker  and  Kays  [l2)  for  0  and  by  Lai  (juQ  for  0. 

The  eigenvalue  problem, 


F(  Yj  ,yU,a,M,T)  -  0 

defined  by  Iqs.  (27)  and  the  boundary  conditions  (26)  at  x  -  ♦£,  can  be  solved 
by  the  Oalerkin  wthod  described  previously.  The  expension  functions  for  6  are 
the  as  those  for  v.  Computations  have  been  carried  out  using  two  and  three 

terms  in  the  series  (9x9  determinant)  for  Yj  ->  1, -  Oj  and  Yj  -  $,  JK  -  0, 
and  0.2  for  various  values  of  N,  and  for  Yj  -  ft-  \  for  various  values  of 
Ray*  Hie  results  are  tabulated  in  Table  3  and  are  shown  graphically  in  Fig.  5. 

For  the  isothermal  case  M  -  0  the  results  are  within  U<  of  the  re¬ 
sults  of  Chandrasekhar  [2,3]  J  and  for  N  /  0  but  yj^-1  they  are  within  1*  of 
those  given  by  Lai  [ll]  .  Further,  the  mtacimm  change  in  the  critical  Taylor 
yjnimhffr  between  the  two  and  three  term  approximations  is  of  the  order  of 
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The  critical  values  of  a  given  in  Table  3  have  only  been  determined  up  to  +O.Ob. 

CONCLUSIONS 

The  present  results  indicate  that  the  Qalerkin  Method  in  conjunction 

with  simple  polynomial  expansion  functions  can  yield  accurate  results  for  a 

variety  of  stability  problems.  Further,  once  the  integrals  corresponding  to 

2  M 

such  operators  as  (DD*  -  a2)  u,  (DD*  -  a2)  v,  and  (D*D  -  a2)  0  are  evaluated 
a  number  of  different  stability  problems  can  be  treated  with  only  minor  addi¬ 
tional  computations. 

In  Fig.  1  the  variation  of  T  yj 2  with  Yj  for  yW.  -  0  is  shown.  Also 
the  theoretical  results  of  Chandrasekhar  and  Elbert  [5]  and  Kirch gas sner  [73 
as  well  as  the  experimental  results  of  Taylor  [l],  Lewis  [2l*],  Donnelly  [25>3 
and  Caldwell  and  Donnelly  [26]  are  indicated.  In  general  the  agreement  is 
excellent.  The  hysteresis  effect  indicated  by  Lewis  appears  to  be  incorrect. 

See  Caldwell  and  Donnelly.  It  is  interesting  to  note  that  Taylor's  correction 
formula  for  gap  size  is  indeed  very  accurate  for  Yj  >  0.8  but  does  begin  to 
deviate  fairly  rapidly  from  the  present  results  for  Yj  <  0.8.  From  a  practical 
point  of  view  the  present  results  indicate  the  correction  in  T  in  comparing 
small- gap  theoretical  (  VJ  ■  1)  values  with  small-gap  experimental  ( Yj  near  1) 
results.  Thus  for  ■  1,  JJ~.m  0,  T  ■  3391  and  for  Yj  ■  .95  (a  reasonable 
value  for  small-gap  experiments),  ■  0,  T  ■  3511  which  is  a  change  of  3.5^. 

As  mentioned  earlier  for  yii-  0,  Trj2  -  (2  )2(l+yj2)  approaches  a  value 

of  approximately  5 00*  as  Vj-eO. 

In  Fig.  2  the  variation  of  Ty'j2  with  Yj  for  *  0.2,  0,  -0.125,  and 
-0.25  is  shown.  Finally  the  regions  of  stability  in  a  (/i^R^d/ y)  )  (d/R^)^  vs 
( n^d/^  )(d/R1)i  plane  are  depicted  in  Fig.  3  for  -  1.0,  0.7,  and  0.5. 
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It  is  well  known  for  the  small -gap  problem,  Yj  -*  1,  that  eren  if 
S~l(r)  is  replaced  by  its  average  value  accurate  results  are  obtained  for 
JU.^0.  Lin  [27]  has  suggested  that  it  should  also  be  possible  to  do  this 
for  the  finite  gap  problem.  However,  some  care  must  be  taken  in  choosing  a 
suitable  average.  It  is  natural  to  consider  the  velocity  rather  than  the 
angular  velocity,  and  to  replace  V(r)  by  its  average  over  the  gap,  V.V  To 
obtain  a  suitable  average  angular  velocity,  Vav  must  be  divided  by  a  suitable 
length  scale  independent  ot ^/U. .  The  requirement  that  the  average  angular 
velocity  be  the  same  as  the  angular  velocities  of  the  cylinders  in  the  limit 
JA  1,  leads  to  the  choice  of  RQ  for  this  length.  Thus  g(x)  in  Eq.  (10)  is 

replaced  by  Ax  +  ^/(l^-^)  where  Yj0  -  (l+>j)/2  and  -  (l-^)/(-lnh  )j 
and  we  can  anticipate  that 

TUj  ♦  — - )  ■  Tgy  (29) 


should  depend  only  upon  Yj  for  ^  }  0.  Actual  computations  show  that  this 
dimensionless  parameter  gives  the  correct  value  of  T  within  2%  for  a  range  of 
Yj  from  0.2  to  1.0  and  for  ^  0.  In  addition,  it  was  found  that  the  critical 
value  of  TftV  is  relatively  insensitive  to  variations  in  YJ  .  By  using  its  value 
at  Yj  ■  1,  namely  1706,  we  obtain  from  Eq.  (29)  after  substituting  for  A^  and 
from  lqs.  (12)  the  extremely  simple  formula  for  the  critical  value  of  T 

T  - _ 1708(1-  >j2) _  (30) 

which  is  correct  within  1^  in  the  range  0.5  ^  j  $  JA}  0.* 

*After  this  work  was  completed  the  authors  learned  that  Rintel  and  Lieber  [33} 
have  also  suggested  an  averaging  process  for  computing  the  critical  Taylor 
number  for  the  finite  gap.  Their  averaging  is  considerable  more  complicated 
and  does  not  appear  to  yield  any  more  accurate  results. 
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For  the  Dean  stability  problem  the  variation  of  (VBd/ i)  )(d/Rg)^  with 
Yj  is  shown  in  Fig.  U.  Of  particular  interest  mathematically  is  the  non¬ 
monotonic  behavior  of  P  with  an  increasing  number  of  terms  in  the  series  ex¬ 
pansion.  For  example  for  Vj  -  1  and  a  *  3*96  the  result  36.86  for  M  ■  2  is 
closer  to  what  we  might  consider  the  "exact"  value  35.9b  than  the  result  37.26 
for  M  -  3.  At  first  sight  these  results  seated  surprising  end  were  consequently 
checked  by  hand  computation.  However,  if  we  bear  in  mind  that  these  are  nonself- 
adjoint  eigenvalue  problems  we  realize  that  we  cannot  expect  necessarily  mono- 
tonic  behavior  of  P  with  increasing  M.  Indeed,  this  result  shows  that  the  well 
known  result  for  self-adjoint,  positive  definite  eigenvalue  problems  that  the 
approximation  to  the  eigenvalue  is  improved  with  increasing  M  does  not  hold  in 
general  for  nonself -ad joint  eigenvalue  problems. 

Finally,  we  consider  the  effect  of  a  radial  tenperature  gradient  on 
the  stability  of  Couette  flow.  As  shown  in  Fig.  5,  positive  and  negative 

temperature  gradients  are  destabilizing  and  stabilizing  respectively.  For yU. 

2 

not  near  yj  it  is  difficult  in  practice  to  achieve  values  of  N,  without  vio¬ 
lating  the  assumptions  we  have  made,  large  enough  to  change  significantly  Ihe 
critical  Taylor  number  from  that  predicted  by  isothermal  theory. 

However,  A^  will  be  small  for  JX.  near  /j  and  hence  it  will  be  pos¬ 
sible  to  have  large  values  of  N.  In  this  case  the  Rayleigh  number  plays  an 
important  role.  A  graph  of  the  critical  Taylor  number  vs.  Ray  number  for  * 
ytO-*  Yj  (N  ■  i)  is  given  in  Fig.  6.  As  can  be  seen  from  this  figure,  for 
sufficiently  large  Ray  numbers  it  is  possible  for  the  flow  to  be  unstable  even 
though  T  is  negative,  that  is  >  Yj  ^ .  This  was  first  pointed  out  by  Tih  [loj 
and  Becker  and  Kaye  [l2]  for  the  small-gap  problem.  Further,  Yih  observed  for 
the  small-gap  problem  as  is  the  case  here,  that  von  Kamaa 

-  18  - 


's  condition  based 


on  an  inviscid  fluid  (that  the  product  of  the  deceit?  |  and  the  square  of  the 
circulation  P  increasing  outward  is  a  necessary  and  sufficient  condition  for 
stability)  no  longer  holds  when  viscosity  is  taken  into  account.  This  can  be 
seen  from  the  plot  of  T  vs.  N  for  (n  -  i).  Pig.  7.  The  points  in 

the  lower  left  hand  quadrant  correspond  to  ^  decreasing  and  P  increasing 
radially.  Thus  for  suitable  Prandtl  numbers  one  can  find  flows  which  are  un¬ 
stable  even  though  ^  p  2  increases  outward  and  the  flow  is  stable  according  to 
the  inviscid  theory. 

In  a  manner  similar  to  that  used  in  deriving  Bq.  (30),  and  asking  use 
of  the  rather  insensitive  dependence  of  the  differential  operators  on  Vj  for 
J  4  1  it  is  possible  to  derive  an  "average-type"  stability  criterion  for 
the  thernal  problaa.  For  yU  >  0,  £  £  ^1  this  relation  is  given  by 

Tav  ♦  (Bsy)^  -  1708  (31) 

where  T^  is  defined  by  Bq.  (29)  and 

(»ay)aT  -  8ay(*1  ♦  .  -i— >2  •  (32) 

This  result  reduces  to  that  given  by  Becker  and  laye  [12]  for  h-^1. 
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Table  1.  Critical  Taylor  matbero  and  corresponding  values  of  a  for  various 
assigned  values  of  Jl  -  n2/  D-i  and  Y)  ■  E^/Rg.  Subscripts  denote 
the  ranber  of  tens  used  in  the  apprmriWlng  series.  The  results 
given  in  eoluan  C  are  those  given  by  fThniilreseMm  [2}  and 
Chandrasekhar  and  Albert  [5J. 


z1 

a 

(i  f)2 

af)3 

-a 

N 

CM 

fH 

C 

V*3  TA  V** 

1.00 

3.12 

1750.0 

1708.5 

1708.1* 

1708.1 

1.021* 

1.000 

1.021* 

o.5 

3.12 

2331.U 

2276.1* 

2276.1 

2275.3 

1.021* 

1.000 

1.021* 

0.25 

3.12 

2792.7 

2727.1 

2725.9 

2725.3 

1.021* 

1.000 

1.025 

0 

3.12 

31*75.0 

3391*.  6 

3391.3 

3390.3 

1.023 

1.001 

1.025 

-0.125 

3.13 

3952.9 

3863.0 

3856.6 

— 

1.023 

1.012 

1.025 

-0.25 

3.13 

1*575.8 

1*1*71*. 1* 

1*1*63.1* 

1*1*62.6 

1.023 

1.003 

1.025 

-o.5o 

3.20 

6589.1* 

61*62.7 

61*18.1* 

61*17.1 

1.020 

1.007 

1.026 

-0.70 

3.3U 

9718.1 

9598.3 

91*35.6 

91*33.0 

1.012 

1.017 

1.030 

-0.90 

3.70 

15537 

15559 

11*952 

11*91*0 

0.999 

1.01*2 

1.039 

-1.00 

l*.oo 

1961*1* 

1971*5 

18703 

18680 

0.995 

1.055 

1.050 

0.9025 

3.12 

1663.5 

1621*.  2 

1621*.  2 

1.021* 

1.000 

1.021* 

0.70 

3.12 

1869.3 

1825.3 

1825.1 

1.021* 

1.000 

1.021* 

0.50 

3.12 

2128.7 

2078.7 

2078.1* 

1.021* 

1,000 

1.021* 

0.20 

3.12 

268U.9 

2622.9 

26a. 6 

1.021* 

1.001 

1.021* 

0 

3.13 

321*5.  i* 

3172.1* 

3168.8 

1.023 

1.001 

1.021* 

-0.125 

3.11* 

3727.2 

361*5.5 

3638.8 

1.022 

1.002 

1.021* 

-0.25 

3.15 

U368.1 

1*276.1 

1*263.1* 

1.021 

1.003 

1.025 
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Table  1.  (Continued) 


n 

a 

<»)!>2 

(t^2)3 

(T*]2)u 

c  t2/t3  Ty^  t2Ai* 

0.90 

0.81 

3.12 

1579.0 

151*1.9 

151*1.9 

1.021* 

1.000 

1.021* 

o.5o 

3.12 

1931*.  9 

1890.1 

1889.6 

1.021* 

1.000 

1.021* 

0.20 

3.12 

21*70.9 

21*15.2 

21*13.7 

1.023 

1.001 

1.021* 

0 

3.13 

3021*.  7 

2958.9 

2955.0 

1.022 

1.001 

1.021* 

-0.125 

3.H* 

3511.7 

31*38.1 

31*30.8 

1.021 

1.002 

1.021* 

-0.25 

3.15 

1*175.0 

1*092.8 

1*078.1* 

1.020 

l.OOl* 

1.021* 

0.80 

0.61* 

3.12 

H*l5.7 

1383.6 

1383.3 

1.023 

1.000 

1.023 

O.liO 

3.12 

1710.8 

1673.9 

1672.1* 

1.022 

1.001 

1.023 

0.20 

3.13 

J2068.1* 

2021*.  6 

2022.8 

1.021 

1.001 

1.023 

0 

3.13 

*2609.6 

2557.8 

2553.U 

1.020 

1.002 

1.022 

-0.125 

3.16 

3112.8 

3055.3 

301*6.8 

1.019 

1.003 

1.022 

-0.25 

3.17 

381*0.5 

3778.1 

3759.9 

1.016 

T.005 

1.022 

0.70 

0.1*9 

3.12 

1260.2 

1233.5 

1233.0 

1.022 

1.000 

1.022 

0.35 

3.13 

11*1*0.2 

11*10.6 

11*09.7 

1.021 

1.000 

1.022 

0.20 

3.13 

1699.2 

1665.8 

1661*.  2 

1.020 

1.001 

1.021 

0 

3.1U 

2229.5 

2190.1* 

2185.9 

1.018 

1.002 

1.020 

-0.125 

3.15 

2760.2 

2718.7 

2708.5 

1.016 

l.OOl* 

1.019 

-0.25 

3.19 

3597.0 

3556.7 

3532.1 

1.012 

1.007 

1.018 

-o.5o 

3.56 

7821.9 

7586.6 

1.031 

— 
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Table  1.  (Continued) 


7 

y* 

a 

CM 

CM 

ST— 

IH 

»f>3 

C 

V*3  TyTjj  VT4 

0.60 

0.36 

3.13 

1112.9 

1091.5 

1091.0 

1.020 

1.001 

1.020 

0.2$ 

3.13 

1272.9 

1249.7 

12lt8.6 

1.019 

1.001 

1.019 

0.10 

3.  lit 

1581.5 

1555.8 

1553.2 

1.017 

1.002 

1.018 

0 

3.15 

1883.5 

1856.6 

1851.5 

1.015 

1.003 

1.017 

-0.125 

3.17 

21*66.9 

21*36.6 

2it2it.8 

1.015 

1.005 

1.017 

-0.25 

3.2lt 

3512.6 

3188. 2 

3453.2 

1.007 

1.010 

1.017 

o.5o 

0.25 

3.1JU 

973.88 

958.35 

957.22 

957.8 

1.017 

1.001 

1.017 

0.20 

3.1U 

I05it.it 

1038.5 

1037.0 

1.016 

1.001 

1.017 

0.10 

3.15 

1263.2 

12it6.2 

1243.7 

1.015 

1.002 

1.016 

0 

3.16 

l57lt.7 

1554.5 

1549.6 

1550. 

1.013 

1.003 

1.016 

-0.125 

3.19 

2250.lt 

2230.3 

2216.4 

2219 

1.009 

1.006 

1.015 

-0.25 

3.32 

377U.5 

3752.0 

3693.9 

3701 

1.006 

1.015 

1.022 

-0.50 

It.  78 

— 

13553 

13323 

13022 

— 

1.017 

— 

o.Uo 

0.16 

3.15 

81t5.7lt 

833.21 

831.81 

1.015 

1.002 

1.017 

0.08 

3.16 

1025.6 

1011.5 

1009.1 

1.014 

1.002 

1.016 

0 

3.17 

1298.6 

1283.5 

1279.1 

1.012 

1.003 

1.015 

-0.125 

3.21* 

2253.0 

2164.2 

2146.5 

i.ola 

1.008 

1.049 

-0.25 

3.76 

5179.1 

U9U2.5 

4830.6 

1.048 

1.023 

1.072 

0.30 

0.09 

3.19 

728.00 

716.56 

714.81 

- 

1.016 

1.002 

1.018 

0.05 

3.20 

846.17 

833.10 

830.70 

1.016 

1.003 

1.039 
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Table  1.  (Continued) 


V 

A 

a 

(t  rj2)2 

(T>j2)3 

(t  rj\ 

IgA  3  Vt4  V*4 

0.30 

0 

3.20 

10$8.7 

1043.1 

1040.0 

1.01$ 

1.004 

1.018 

-0.12$ 

3.44 

2623.8 

2$25.4 

249$.  2 

1.039 

1.012 

l.o$l 

0.20 

0.0b 

3.22 

624.36 

609.12 

606.72 

1.02$ 

1.004 

1.029 

0 

3.23 

849.24 

833.  ItO 

829.14 

1.019 

1.00$ 

1.024 

0.10 

0.01 

3.29 

$4i.ii 

$13.87 

$09.10 

1.0$3 

1.009 

1.062 

0 

3.30 

673.59 

6$7.09 

650.00 

1.02$ 

1.001 

1.036 

to*  wh'  j  pi  inp*  i 


Table  2.  Critical  values  of  R(d/fcg)*  -  (V-d/*>  Xd/Rj.)*  aw)  oorraepoodia*  values 
of  a  for  various  assigned  values  of  W  .  Subscripts  daaota  tfca  aabv 
of  Uraa  used  in  the  apprcadjaating  series. 


1 

a 

[iXd/Rg^Jg 

[iKd/Rg)*^ 

R2/R3 

VRi* 

— f — — — 

1.0 

3.96 

36.81* 

37.31 

35.90 

0.989 

1.039 

1.025 

0.95 

1*.02 

37.19 

37.70 

36.26 

0.986 

1.01*0 

1.025 

0.90 

l*.o6 

37.72 

38.30 

36.79 

0.985 

1.01*1 

1.021* 

0.80 

l*.l6 

38.82 

39.51 

37.88 

0.983 

1.0l*3 

1.021* 

0.70 

U.  21* 

1*0.11 

1*0.96 

39.19 

0.979 

l.Ql*5 

1.023 

0.60 

U.32 

1*1.61* 

1*2.73 

1*0.79 

0.975 

1.01*6 

1.021 

0.50 

1*.  1*1 

1*3. 1*7 

1*1*.  91 

1*2.79 

0.968 

1.050 

1.017 

0.2*0 

2*.  1*6 

1*5.65 

1*7.70 

1*5.18 

0.957 

1.053 

1.011 

0.30 

l*.5l 

1*8. 2*7 

5i.l*7 

1*8.75 

0.91*2 

1.056 

0.995 

0.20 

U.57 

51.72 

57.00 

53.60 

0.908 

1.063 

0.961* 

0.10 

1*.61* 

57.03 

65.91 

60.92 

0.86$ 

1.062 

0.936 

26 
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Table  3*  Critical  Taylor  numbers  and  corresponding  values  of  a  for  various 

D  assigned  values  of  ^  »yU  >  and  I  or  Say.  Subscripts  d suets  the 

number  of  terms  used  in  the  approximating  series. 


0  1 

A 

a 

N 

Ray 

T2 

T3 

n  1-° 

0 

3.12 

1.0 

2223. U 

2177.3 

1.0212 

Li 

o.5 

3712.1 

2653.3 

1.3990 

n 

0 

3U7U.7 

339U.6 

1.0235 

U 

-o.5 

U829.U 

U705.9 

1.0262 

[i 

-o.75 

5992.0 

5826.3 

1.028U 

0 

-1.0 

787U.9 

7631.9 

1.0318 

fi  °-s 

0 

3.15 

1.0 

U659.9 

U609.O 

1.0110 

Li 

o.5 

5356.9 

5295.3 

1.0116 

c 

0 

6296. U 

6218.2 

1.0126 

-0.5 

7628.7 

7522.2 

1.01U2 

1 J 

-0.75 

8526.3 

8397.3 

1.015U 

0 

-1.0 

9656.6 

9U9U.8 

1.0170 

y  °.5 

0.2 

3.15 

1.0 

2895.2 

2852.3 

1.0150 

0.5 

3U35.U 

3382.8 

1.0155 

D 

0 

U221.9 

U15U.2 

1.0163 

n 

-o.5 

5U71.3 

5376.8 

1.0176 

Ll 

-0.75 

6i*i7.i 

6300.0 

1.0186 

\ 

-1.0 

7750.7 

7598.0 

1.0201 

li 

I: 

II 


Table  3.  (Continued) 

)  f.  . 

N 

Ray 

T2 

t3 

V*3 

0.3  0.23  3.13 

16000 

-3996.7 

-UoUl.2 

.989 

12000 

-2010.2 

-2038.1 

.977 

7828 

0 

Uooo 

19375 

1879.2 

1.031 

0 

3898.3 

3833.6 

1.017 

-Uooo 

3850.6 

5778.5 

1.012 

-7000 

7309.1 

7231 J. 

1.011 

-60  -40  -20  0  20 
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V&«.  3,  H»  rtobl*  sad  aaatotto voglaaa  la  a  (/l^d/* 
(AA</i>  X^)*  plaaa  f*  wrUaa  tolaoa  of  7. 


APPENDIX 


All  of  the  integrals  that  are  needed  to  use  the  Qalsrkln  Method  with 
the  polynomial  functions 

%(*)  -  (x2  -  *)2  x“-l  ,  T„(x)  -  e^x)  -  (x2  -  *)  (Al) 

can  be  expressed  as  linear  combinations  of  the  following  three  functions: 


(A3) 


ri 

P(m,n,  )  -  ^(x2  -  ln^  dx 

A 

(-£/nJk 

-  S(mfn)ln q  +  £,S(n*l,n)  ♦  ^  Z-  — jy -  S(M*,n)  (AU) 
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Hare  a  and  n  are  zero  or  positive  integara,  p  is  a  positive  integer,  and 

•  (is) 

The  series  converge  for  0  <  >j  <  and  converge  for  Yj  -  0  whenever  the  inte¬ 
grals  do.  If  Vj  -  *1,  then  £  -  0,  and  the  series  either  terminate  after  one  tern 
or  are  identically  zero.  The  series  converge  rapidly  for  Yj  in  the  range  of 
interest.  For  example,  to  evaluate  F(l,l,l,i)  and  F(l,l,l,0.l)  correct  to  six 
significant  figures  requires  8  and  27  terms  respectively. 

Let 


L  -  DD*  - 


d  /  d  + 

SE(S 


L*  -  D*D 


then 


(’’m^n)  -  S(ra*n-2,3) 
(u  x,v  x)  -  -  l/lljO 
(u1,v2)  -  0 
(upVj)  -  -  l/$0h0 
(u2,v2)  -  -  l/$0h0 
(u2,V3)  -  0 
(uyu^  -  -  1/73920 
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<3  vb>v  - 

VJQ S(»*n-2,2)  ♦  £S(»m-l,2) 

()  *L»*1>  - 

V* 

(3V2} " 

E/BUo 

(M  * 

V81*0 

* 

£-Aoo8o 

(^y2,v2)  - 

y*  1*0 

(T  v2»v3^  " 

<£A0080 

( 1  ▼2»V  * 

^oAooeo 

(}v3,v3)  - 

*]oAoo8o 

0V3»VU)  “ 

<6/88701 

(JTu*V  “ 

y8870U 

- 

^0S(m+n-2,U)  ♦  E  S(a»n-l,U) 

- 

y630 

(^ui,^)  - 

6/27720 

(^«l#u3)  “ 

^27720 

( )  uliuU)  ■ 

£,/li80U80 

0  Vtt2}  - 

^27720 

( }  u2»u3^  " 

SA80U80 

( J  u2»ulP  “ 

yU80ltf0 

(}u3,u3)  - 

yii80l;80 

0U3»V  ■ 

6/5765760 

<5vV ' 

^5765760 

Ov»*v  - 

l^S(m+n-2,3)  ♦  £s(»*n-l,3) 

Ovi»v " 

-yu* 

(^▼l»«2)  - 

-£/5oUo 

(}rl»u3)  “ 

-rje/soho 

Ovi»V  ■ 

-5/75920 

(}T2»U2>  - 

-  'Jo/^OliO 

()  t2>«3)  - 

-  5/73920 

(}v2»V  " 

-  ^73920 

<ir3>*3>  " 

-  >]o/73920 

Gv3»V  " 

5/768768 

<}Tu»V  - 

-  ^768768 

C  j  -.i»vn) 

^  os(®fn”2»3)  ♦  2  ^0£S(**-n-l,3)  ♦  £„2S(*m,3) 
-  l/U*0  -  £2/$oUo 

( 3  2*i,vi) 

(3  ^1^2) 

-  d  ^2520 

(  )\,*3) 

-^2/50U0  -  £i/73920 

('j2U2,T2) 

-  ^|/50l*0  -  5.2/73920 

(}2u2,t3) 

-  5^36960 

(  ^  2U3,y3) 

-^|/73920  -  62/768768 

(  3"\»un> 

-  F(**n-2,3,1,  rj ) 

(  3"2%>Tn) 

“  F(»*n-2,3,2,  rj ) 

J  )  “  P(»*n-2,3,>J) 

-  la.  - 


I 

0 

[; 

L 

L 

L 

L 

D 

C 

0 

D 

0 

L 

D 

L: 


L 


(}Tl» Lvi)  -  “  V3  "  e2^0*2*1**)) 

(jv^Lvg)  -  -  d/30  -  £2F(l,2,l,  Vj  ) 

(  J  t1»i,v3^  "  ”  %/6°  -  £2F(2,2,1,  *j ) 

(}vl»LV  "  "  ^Z280  -  e2F(3,2,l,  YJ ) 

(|v2,Lv2)  -  -  rjj 20  -  £2F(2,2,1 ,»J) 

<7y2»Lv3>  “  -  -  £2F(3,2,1  ,/j) 

(^▼2»LV  "  "  3  7o/560  -  £2F(U,2,1,>)) 

(]v3,Lv3)  -  -  11  ^1680  -  £2F(U,2,1,  >)  ) 

(}v3,Lvu)  -  -  ll£A0080  -  £^F(5,2,i,  Vj) 

(J  -  -  23  ^20160  -  £?F(6,2,1  ,rj) 

(  |  u^Li^)  ■=  -  2  ^qA05  -  e2Ho,k,i,fj) 

(  Ju^ug)  -  -  e/l260  -  £2f(1,U,1  ,YJ) 

(  ;ui,lu3)  -  0  -  e2 ?(2,h,ity 

(  }  U;l,Luu)  -  -  e  /mo  -  £2F(3,h,l>  /j  ) 

(ju2,Lu2)«-  ^q/630  -  £2F(2,U,1^) 

(]u2,Lu3)  -  -  e/11088  -  £.2?0,k,i>rj) 

(}u2,Luu)«-  ^q/13860  -  e^CU.UA,^) 

(  }u3,Lu3)  -  -  >)o/92U0  -  £?F(U,U»1»  ^  ) 

(  }  u3,Lu^)  -  -  31  e/2882880  -  £,2F(S,1;,1i  A  ) 
(Ju^Luj^)  -  -  17^oAW4lU*0  -  £?F(6,U,l*jJ> 

( jTl»L^l)  "  -  >]o/3 

(  ^v1,L*v2)  -  -  6/30 
( j  rltL*r3)  -  -  ^(j/60 

-  1*2  - 


(}t2,L*t2)  -  -  Yjo/20 
( )  t2,L#v3)  -  -  6/8I1 

(  }  t3»1#t3)  “  -  11  ^  q/1680 

( }  -  U  rjQ/$  -  3  C2  [l2P(2,2,l,  *j )  -  *(0,2,1,  >1 )] 

(  -  3t/3$  -  3£2  [lfip(3,2,l,)j  )  -  jf  (1,2,1,  >) )] 

(  }  Vl^)  •  -  3  £2  [21F(U,2,1,>|  )  -(7/2)P(2,2,l,  )  ♦(lA6)P<0,2,l,>|  )] 

-  dAUO  -  3£  2  [27F(5,2,1,*|>  -0L1/2)P(3,2,1,J|)  ♦  (3A6)P(l,2,l,  h  )] 

(  }  U2,l2u2)  -  y7  -  3  £2  [2QP(U,2,1*  *]  )  -  3F(2,2,1,  >|  )] 

(  ^U2^2u3)  -  £/60  -  3£2  [25P(5,2,l,y| )  “(?/2)P(3,2,l,Vj )  ♦  (Vl6) P(l,2,l,  KJ  )] 

(  }  iig,!2^)  -  ^0/8U  -  3  L2  [31P(6,2,1,>] )  -(13/2)f(U,2,1,>|)  +  (}A6)F(2,2,l,ir|  )] 

(  Ju3,l2u3)  -  3  tj'JiUo  -  3  e2  Q3»(6,2,lf  >1 )  -  6F(U»2,1,V|  )  ♦  (l/8)p(2,2,l,  >J  )J 

0  u3,l2^)  -  59  £ A8U8O  -3  e2  [36P(3,U,1,  >)  )  ♦  10F(3,3,1,  >| )  ♦  (l/2)F(3,2,l,  Vj )  ] 

(ju^lA^)  -  13  ^  o/3696  -  3t2  [U2F(U,U,1,  >j)  ♦  UP(U,3A,»)  )  ♦(l/2)F(U,2,l,  >j)] 

It  1b  possible  of  course  to  derive  fonulaa  for  the  pneNl  lamer 
products.  For  example 

(  J v*,Lrn)  -  jm(«*-l)S(«*n-2,l)  -  (i)(»-l)(m-2)S(m*n-U,l)J 

.♦  d  {  (■H)2S(*m-l,l)  -  (*)(o-l)%(am-3A)j 
-  £2F(*m-2,2.l»>l) 


